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We discuss flat-band surface states on tfie (111) surface in the tight-binding model with nearest- 
neighbor hopping on the diamond lattice, in analogy to the flat-band edge states in graphene with 
a zigzag edge. The bulk band is gapless, and the gap closes along a loop. The verge of the flat-band 
surface states is identical with this gap-closing loop projected onto the surface Brillouin zone. When 
anisotropics in the hopping integrals increase, the bulk gapless points move and the distribution of 
the flat-band states expands in the Brillouin zone. Then when the anisotropy is sufficiently large, 
the surface flat bands cover the whole Brillouin zone. Because of the completely flat bands, we can 
construct surface- state wavefunctions which is localized in all the three directions. 
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I. INTRODUCTION 

Flat bands have been studied particularly in the con- 
text of possible ferromagnetism driven by interactions, 
as was proposed by Lieb^, and successively Mielke and 
Tasaki^^. On the other hand, from the research on 
graphene^ it is known that the tight-binding model with 
nearest-neighbor hopping on a honeycomb lattice with a 
zigzag edge exhibits flat-band edge states^, and its ori- 
gin is topologically interpreted^. In the dispersion of a 
graphene ribbon with zigzag edges, the flat bands appear 
between the wavenumbers corresponding to the projec- 
tion of Dirac points at K and K'. In contrast, there is 
no fiat band in the edge states in the graphene ribbon 
with armchair edges because in the projection of the dis- 
persion, Dirac cones at K and K' overlap each other. 
In three-dimensions, the diamond lattice, i.e. a three- 
dimensional analogue of honeycomb lattice, exhibits flat- 
band surface states, as theoretically found in the nearest- 
neighbor tight-binding model^. 

When the hopping of the tight-binding model on the 
honeycomb lattice becomes anisotropic, the Dirac points 
in the bulk BZ move away from K and K'. Moreover when 
the tight-binding model has sufficiently large anisotropic 
hopping integrals, the two Dirac points meet and the bulk 
dispersion relation is linear in one direction and quadratic 
in the other—. In that case, the flat-band edge states 
cover the whole one-dimensional (ID) BZ^^. With the 
increase of the anisotropy, the bulk becomes gapped while 
the complete flat band remains tn the edge BZ. 

In this paper we focus on surface flat bands in the tight- 
binding models on the diamond lattices with (111) sur- 
face. When the model has no anisotropy, the gap closes 
along a loop in the bulk BZ. If the Fermi energy is set to 
be zero, which corresponds to the case with particle-hole 
symmetry, the bulk Fermi surface (FS) coincides with 
this loop. An anisotropy in the nearest-neighbor hopping 



integrals deforms the loop. Similarly to the honeycomb- 
lattice model, the surface flat bands are formed in the 
k region surrounded by the projection of the FS loop. 
When the anisotropy is sufficiently large, the FS loop 
vanishes, and consequently the flat-band surface states 
cover the whole 2D BZ. 

Furthermore, because the edge/surface band is com- 
pletely flat over the entire BZ in the models on honey- 
comb/diamond lattices, any linear combinations of the 
edge/surface states remains eigenstates. Thereby we can 
construct the edge/surface states which are spatially lo- 
calized in all directions. In the diamond-lattice model 
we also find a topological transition of the bulk FS by 
changing anisotropy; the loop of the bulk FS changes its 
topology by varying the hopping integrals. 

The organization of the paper is as follows. In Sec. II 
we review how the flat-band edge states of the tight- 
binding model on the honeycomb lattice evolve as we 
change an anisotropy. We discuss analogous behaviors 
of surface states of the model on the diamond lattice in 
Sec. III. In Sec. IV, we show how the behaviors of the 
edge/surface states shown so far are explained by topo- 
logical argument. Section V is devoted to a calculation 
of the states which is localized along the boundary, for 
both the honeycomb- and diamond-lattice models. We 
summarize our results in Sec. VI. 



II. HONEYCOMB LATTICE 

We flrst review the flat-band edge states on the 
honeycomb-lattice structure shown in Fig. [IJa), and 
study the completely flat band, which has been studied 
in Ref. 11. We consider a tight-binding Hamiltonian on 



this lattice, 
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where the suffix h represents the honeycomb lattice, tij 
is the hopping integral along the nearest-neighbor bond 
vector Ta, and q (q) is the annihilation (creation) oper- 
ator of the electron. We treat the hopping integrals Uj as 
real positive parameters, and they are labeled with the 
vectors Ta as ta- The bulk Hamiltonian matrix H\^\^{k) 
at wavevector k is given as 



i^hb(k) 







Eli^^e^^-i 



Eti^^e-^^-i 







(2) 



As we see in the following, for t > 2 fiat-band edge states 
on the zigzag or Klein edges completely cover the BZii. 
It occurs when the bond with hopping t is perpendicu- 
lar to the edge for the zigzag edges, or when the bond 
with hopping t is not perpendicular to the edge for the 
Klein edges. We calculate dispersions in Fig. [1] (c) for 
ti = 1 and ti — 2.2, ^2 = fe = 1. The flat-band 
edge states are separated completely from the bulk for 
t\ = 2.2 (Fig. [TJc)). To explain this behavior, we solve 
the Schrodinger equation in the semi-infinite geometry 
with a zigzag boundary. The zigzag edge is assumed to 
be perpendicular to the bonds with hopping integral ti. 
We express the wavefunction |^(A:)) as 



\^{k)) =J2i^,{k)\Mk)) +b,{k)\B,{k))), 



(5) 



where the suffix b means the bulk. Ti=i^2,3 are expressed 

as Ti = (0,1), T2 = (-^,-|),T3 = (^,-|), where we 
put the length of the nearest-neighbor bonds as unity. 
For simplicity ti are assumed to be positive. The primi- 
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tive vectors s.i=i^2 are ai = (^, |),a2 = ( 
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). The 



reciprocal primitive vectors are Gi = 27r|(^, |),G2 = 

^^3V 2 ' 2r 

We ffist note that the bulk Hamiltonian i!/^hb has 
particle- hole symmetry: osHi^h^^s = — ^hb- Therefore, 
if \ip) is an eigenstate with the eigenvalue E^ crslip) is an 
eigenstate with energy —E. The eignenvalue is given by 



where i denotes the index for unit cells containing two 
sublattices, A and B, counted from the edge {i = 1), 
k is the wavenumber along the edge, and a^{k) {b^{k)) 
denotes the amplitude of the Bloch wavefunctions \A-{k)) 
{\B-{k))) for the i-th A (B) sublattice. Acting Hi^ onto 
\A-{k)) and \B-{k)), we have 

{B^{k)\H^\Mk))=h, (6) 

{B,_,{k)\H,,\Mk))=h^tse-^\ (7) 

When the eigenvalue is zero, Hh\'^{k)) = 0, we obtain 



Shb(k) = ± 






-iTi-h 



(3) 



aiit2 + he-''') + tti+iti =0, hi = 0. 



(8) 



Hereafter we put parameters ^2 = H = 1, and ti = t, 
where t is a real positive tunable parameter. The bulk 
dispersion is given as 



Thus the condition for existence of the edge states, i.e. 
normalizability of the wavefunction, is given as 



t2 + fee 



U 



<1, 



(9) 



^hb = 
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The bulk gapless points {k^^ky) are given by the equa- 
±^t and sin I A: * = 0. The equations 



and we have the amplitude of the fiat-band states as 
fe + fee"* 



a^{k) = ai 



tions: cos 



VsK 



give two gapless points in the bulk BZ, and they exist 
for t < 2. The gapless points move with the change 
of the anisotropy t. For t = 1, i.e. the tight-binding 
model of graphene, the upper and lower bands touch at 
K (^^, ^) and K' (-^^, ^), and with the increase 
of t the gapless points get closer along the line ky = ^ 
(Fig. mb)). Around each of the two gapless points, the 
dispersion forms a Dirac cone, and Berry phase around 
each gapless point is tt. Because of this tt Berry phase, 
the gapless points do not disappear as we change t{< 2). 
The bulk gapless points move in the direction perpen- 
dicular to the anisotropic hopping integral t. At t = 2 
the gapless points meet and they annihilate each other at 
kx = (Fig.djb))^^. For t > 2, there are no bulk gapless 
points. 

The evolution of the edge states with the change of 
the anisotropy has been studied in several papers!^— . 
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, bi = (10) 
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from normalization. 



where a-^ 

For ti = fe = fe (graphene model), the wavenumber 

which satisfies the condition (Eq. (J9j)) for existence of the 

edge state is given as ^ < /c < ^ , which agrees with the 

27r well-known fiat band in graphene ribbon with a zigzag 



t2+^3e 



< 
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edge^. In addition, by the relation 

the fiat bands cover the whole ID BZ when fe > fe +fe is 
satisfied. This is realized when anisotropy is sufficiently 
large. These results agree with numerical calculations in 
Fig- [11(c). In Fig.[T](c) we also show results for armchair 
edges and for Klein (bearded) edges. For armchair edges 
there are no fiat-band edge states. For Klein edges where 
Ti-bonds (hopping fe) are not perpendicular to the edge, 
there are fiat-band edge states; if fe > 2 the fiat-band 
edge states cover the entire BZ. These results agree with 
the results in Ref. JT,. 
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FIG. 1: (a) Schematic of the honeycomb- lattice structure. 
The dotted hne shows a choice of unit cell with translational 
symmetry along the edge. The arrows show the directions 
along zigzag, armchair, and Klein edges, (b) shows the first 
BZ of the honeycomb lattice. Dots at the corners show the 
gap less points. The bulk bands become gap less at K and K' 
points, when the hopping integral is isotropic. By increasing 
ti from t2 — t2> — 1, the gapless points move away from K 
and K' points as shown by the arrows. The flat-band edge 
states expand in the BZ as the bulk gapless points move, (c) 
shows the dispersions for ti = 1 and 2.2 at t2 = ts = 1. The 
dispersions are shown for zigzag, armchair and Klein edges. 
In zigzag and Klein edges, flat bands appear at the Fermi 
energy. 



III. DIAMOND LATTICE 

In the previous section, we have seen that the flat band 
of the honeycomb lattice with the large anisotropy covers 
the whole edge BZ. In this section we show surface flat 
bands^ in the tight-binding model of the diamond lattice 
in the same method in the previous model on the honey- 
comb lattice. In the honeycomb lattice, the verge of the 
flat bands is identified with the bulk gapless points. We 
will show the same phenomena for the diamond lattice. 
The Hamiltonian H^ is 



H. 






(11) 



where tij = tji is the hopping integral from site i to j, 
and the suffix d means the diamond lattice. We assume 
that the hopping integrals are real positive parameters. 
The nearest-neightbor bond vectors, rs, are as follows: 



n = i(i,i,i), T2 = |(-1,1,-1), 



T3 



i(-i,-i,i), 



and r4 = |(1,— 1,— 1) (Fig. [2j). The hopping integrals 




FIG. 2: (a) Schematic of a diamond lattice. A (gray) and B 
(black) atoms denote the sublattice. rs represent vectors con- 
necting nearest-neighbor atoms, (b) The BZ of the diamond 



lattice structure within kx,i 

are labeled with Ta as t^ 
i7db(k) is represented as 



>0. 



ta- The bulk Hamiltonian 



i^db(k) 



Et 







E^=i Ue 




-ik-Ti 



(12) 



This also satisfies <j3J^dbC'"3 = —H^h^ and therefore the 
eigenstates have particle-hole symmetry. As is proposed 
in 9, when t^'s are identical, the surface states on the 
(111) surface for this Hamiltonian forms a flat band, 
which partially cover the surface BZ. 

First, we set one of the four hopping integrals as a tun- 
able positive parameter t and the other hopping integrals 
as unity. For t = 1.4, the dispersions for slab geometry 
with (111) surface are shown in Fig.[3l^a)(b). In Fig.[3l^a) 
we set ti = t, t2 = ts = t4 = 1 and in Fig. [3jb) we set 
ts = t, h = t2 = U = I. The left panels of Fig. ^di)(h) 
show the dispersions near zero energy, and the distri- 
butions of zero-energy flat bands in the surface BZ are 
shown in the right panels in (a)(b). Fig. [3fc) shows the 
projections of the bulk FS {E = 0) onto the (111) sur- 
face for (a)(b). By comparing the FSs in Fig. [31(c) with 
the flat-band states (a)(b), we see the correspondence be- 
tween the bulk FS and the distribution of the flat bands 
as is discussed in the previous case of the honeycomb lat- 
tice; the verge of the flat bands corresponds to the bulk 
FS projected onto the surface direction. The bulk FS 
forms a loop in the bulk BZ because Hdh consists only 
of (Jx and (jy. In this case, the FS loop encircles the axis 
parallel to the hopping vector with the hopping integral t. 
Therefore by rotating the surface orientation in (b) from 
(111) to (111), the bulk FS is rotated as well, and corre- 
spondingly we have the same flat-band states as shown 
in (a). 



The bulk dispersion for ti = t, ^2 = ta = ^4 = 1 is 
given by 

2 



E, 



db 



t + COS + COS — ^ h COS 



gj^ ^l+fc. ^ gj^ fcj^ ^ gj^ kz + k^ ,^^3^ 



Therefore, the gapless points £^5, 
equations, t + cos 






■ cos ■ 






sm T +sin ^ 



are given by two 

+ cos ^^±^ = and 

2 I ^^^x 2 — ^ +sin ^^^^'^ = which determine the 
bulk FS are shown in Fig.HJ The FS encircles the V-L line 
in the bulk BZ for t > 1. We see that the FS is getting 
smaller with the increase of ti and shrinks to the F point 
at t = 3, while the flat-band surface states expand with 
the increase of ti. For t > 3 the loop vanishes, and the 
surface fiat band covers the whole surface BZ. 




FIG. 3: Band structure for the diamond lattice model with 
the (111) surface, for 12 — 1^ — 1 and (a) ti — lA and t^ — 1 
or (b) ti = 1 and ts = 1.4. The right panels in (a)(b) show the 
region with the surface flat bands. In (c), the FS {Ehd — 0) 
in the 3D bulk BZ projected on the (111) plane for the cases 
(a) and (b). 
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FIG. 4: The FS [Em = 0) in the 3D bulk BZ projected on 
the (111) plane for t2 = ts = t4 = 1 and ti = 0.6, 1.4, 2.2, 3.0 
are shown. For ti > 1, the loop shrinks with increase in ti, 
and at ti = 3.0 the loop becomes a dot. 



where /ci = k • (t2 — T3) and /c2 = k • (ra — t/^). In the 
Schrodinger equation, i^d|0) can be calculated as 



i^d|0(k)) 



i=2 



iki 



h + U^' 



ik2 



■a^^,h]\B^{k)) 



ts + t4e 



ik2) 



IA(k)) 



Y^ibj-ih + bj{t2e'''' + ig + ^46" 



■ifca^i 



\AAm 



j='2 



From H(i\(j)) = 0, the amplitudes are derived as bi{k) = 0, 
and 



^j+i 



(k) 



t2e 



ts + Ue'' 



U 



-a,.(k). 



(17) 



A. Flat-band states localized at the surface 

Here we show the Bloch wavefunction of the flat-band 
states for the tight-binding model (Eq. (pT]) ) with a (111) 
surface. To this end, we consider the tight-binding model 
on a semi-infinite geometry of the diamond lattice with 
the (111) surface. We derive a wavefunction with zero 
eigenvalue as \<j>) = Ei=i(ai(k)|^i(k)) + 6,(k)|B,(k))), 
where i denotes indices of the unit cell containing A and 
B sublattice sites, counted from the top {i = 1), and 
a^(k)(6^(k)) denotes the amplitude of the Bloch wave- 
function for each sublattice. The matrix element of the 
Hamiltonian i^d is given as 



(B,_i(k)|JTd|A,(k))=ti, 
(B.(k)|ifdlA-(k))=t2e-^'^^ 



ts + t^e 



ik2 



(14) 
(15) 



Namely, the zero-energy states localize near the surface 
for every wavenumber for ii > ^2 + is + ^4, because we 
have 



a^-+i(k) 



aJk) 



tae-^'^i + ^3 + t4e'''^ 



h 



tl 



from Eq. (p!7|) . The second equality in Eq. (p!8]) holds 
only at the F point in the 2D surface BZ. In addi- 
tion, from Eq. (fT8|) a penetration depth A, defined as 
\aj\ ~ e~^l^^ is calculated for each k. Equation (p!8|) 
gives e~A = ^^^ — ^+^^+^6 ^ Hence the penetration 

depth of the surface states is maximum at k = 0, and 
the finite-size effect is the largest. 

From the above discussion, the model can have the 
complete flat band at zero energy when the system is 
semi-infinite. However, if the thickness of the slab is 
finite, there is a small gap around the F point, as shown 



in our numerical results for surface-state dispersions in 
Fig. [5fa)(b). The maximum value of the finite-size gap 
A is estimated as the following. We focus on the F point, 
and assume (^2+^3 + ^4)/^! = a < 1. When the thickness 
of the slab N is large, the surface wave function for the top 
surface \(j)^) and that for the bottom surface \(j)^) can be 
treated separately. The wavefunctions are approximately 



dotted fine, and the straight lines shows the FS at t = 1: 
kx = ±27r or ky = ±27r at kz = 0^ kx = 27r at /c^ = 
and ky = 27r dit kx — 0. The squares at /c^ = and 27r 
are equivalent. At t = 3 the FS shrinks to a point and 
simultaneously the surface flat band extends cover the 
whole BZ. 



given as 



N 



N 



|0*(k)) ~ ^ac"-V.(k)), |/(k)) ~ ^ac^-"|B,(k))(^9) 



n=l 



where c = -{t2e~'^^ + ^3 + Ue'^^)/ti, and a^ = (1 



^)/(l 



,|2Ar 



). Then we obtain the hybridization as 



\]<.)\H\^'0<.))^a'c^h, 



(20) 



which is expected to give the size of the gap due to the 
finite-size effect. At the F point, \c\ becomes maximum 
and therefore the finite-size effect of the energy is largest 
at F, taking its maximum value 



A ~ ti{l - a'^)a 



N 



(21) 



For ti = 3.4, ^2=^3=^4 = 1.0, ha^ - 0.0616 for N = 
20, and ha^ - 0.00504 for N = 40. Our band-structure 
calculation gives A to be A ~ 0.0632 for N = 20, and 
A = 0.005042 for TV = 40, as shown in Fig. [5l Thus 
our estimate for A in Eq. (|2T]) well agrees with numerical 
calculation. 




FIG. 5: Surface dispersion near the Fermi energy at t = 3.4. 
The number N of unit cells to the direction normal to the 
surface is A^ = 20 in (a) and A^ = 40 in (b). The surface BZ 
is the same as that in the left panel of the Fig. [Sja). 



B. Topological transition of the bulk FS 

In this section, a topological transition of the FS in the 
diamond lattice is discussed. We assume that ^2,3,4 are 
fixed to be unity, whereas ti = t is varied. To illustrate 
the change of topology of the FS, we choose the scheme of 
the BZ as shown in Fig. [6l For t < 1 the FS is separated 
in the BZ. The Fermi loop reaches the boundary of the 
BZ, and forms an open orbit. On the other hand, for 
t > 1, the loop is connected in one closed loop in the BZ, 
and the FS shrinks with an increase of t. The transition 
of the topology of the FS occurs at t = 1. The FSs for 
t = 0.8, 1, 1.2 are plotted in Fig. [6] to show the process 
of the transition. In the figure the BZ is shown as the 



{^^^)//i 




t = 0.8 






-2^ kx 
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FIG. 6: The FSs in the 3D bulk BZ for t2 = ts = t4 = 1 
and ti =0.8, 1, 1.2 are shown. At ti =0.8 the FS is separated 
in the BZ. On the other hand at ti = 1.2, the FS becomes 
one closed loop in the BZ. The topological transition occurs 
at ti = 1. At the right side, the FSs around the X point are 
magnified. The change in FS topology occurs at the X point. 



IV. TOPOLOGICAL EXPLANATION FOR 
EXISTENCE OF THE FLAT-BAND STATES 

Here we give a topological reason for existence of 
the fiat-band edge/surface states in the present mod- 
els, based on Ref. 8. For the honeycomb-lattice model, 
for example, we first begin with a bulk system, and 
later we cut the system along one direction by cutting 
the nearest-neighbor bonds, in order to discuss edge 
states. Let y denote the coordinate along which the 
system will be cut. Then, following Ref. 8, we ex- 
pand the bulk Hamiltonian by the Pauli matrices ax^y 
as iJ = hx(k'^^ ky)cFx + hy{k'^^ ky)(jy. Here ky denotes the 
component of the wavevector along the ^/-direction (along 
the edge) , and /c^ is another component of the wavevec- 
tor. If the bulk system is cut along the ?/-axis, k'^ will no 
longer be a good quantum number. Then the criterion in 
Ref. 8 says that if the trajectory of {hx^ hy) for the change 
of k'^ with fixed ky encircles the origin, zero-energy edge 
states exist for the given ky. If not, zero-energy edge 
states will not exist^. 

We apply this criterion to the present models to show 
that the fiat-band boundary states discussed so far are 
fully explained by this theory. For the anisotropic 
honeycomb-lattice models, explanations are given in 
Ref. ill,. For zigzag edges we have 



hx = t cos{ky — k'^) + 1 + cos ky^ 
hy = —tsm{ky — k'^) -\- sin ky. 



(22) 
(23) 



(a) \h 



(b) hy 




FIG. 7: Trajectories of (hx^hy) for the honeycomb-lattice 
model with zigzag (a), with Klein (b) edges, and in the 
diamond- lattice model with the (111) surface (c). In these tra- 
jectories ky, ki and /c2 are fixed while k^ and ks are changed. 



Hence the trajectory is a circle with a radius \t\ centered 
at (1 + cos/c^y, sin/c^y) (Fig. [Tl^a)). The condition that it 
encircles the origin fully reproduce the behavior of the 
flat-band edge states. In particular, for t > 2 the trajec- 
tory encompasses the origin irrespective of the value of 
ky, and existence of the perfectly flat band results. The 
case for the Klein edge is explained similarly, where we 
have 



hx = 1+ cos(ky - k'x) + t cos k'x, 


(24) 


hy = - sin{ky -k'x) + t sin k'^. 


(25) 



with the trajectory shown in Fig. El^b). The flat-band 
edge states extend over the whole BZ when t > 2. 

Moreover, this discussion can be extended to three- 
dimensional model such as our diamond-lattice model. 

hx = I -\- cos ki + cos k2 -\-t cos /cs, (26) 

• ' -• ' ■-• ■ (27) 



hy = sin ki — sin k2 — t sin ks, 



where ki = k • (t2 - T3), /c2 = k • (T3 - T4), and ks = 
k • (t3 — ti). When we cut the bonds parallel to ti, 
ks will no longer be a good quantum number. Then, if 
the trajectory of ihx^hy) by the change of ks encircles 
the origin, there should be zero-energy edge state. For 
t > 3 it holds true irrespective of the values of ki and /c2 , 
and therefore the surface flat band covers the surface BZ 
(Fig. He)). 

This explanation also implies that in these models, the 
verge of the edge/surface states should be identical with 
the projection of bulk FS where the conduction and va- 



lence bands touch each other. Let k* denote a wave num- 
ber where the bulk FS satisfies hx = ^ and hy = 0. In 
our case the honeycomb lattice model has two k* points 
(Fig. [1]). k* in the diamond lattice model forms a loop 
(Fig. [3] and [4j). Their projections onto the edge/surface 
are identical with the verge of the flat-band edge/surface 
states, which is one of the consequences in Ref. 8. Thus 
this topological notion enables us to show existence or ab- 
sence of flat-band boundary states for various boundary 
conditions. When the flat-band edge/surface states exist, 
the region enclosed by k* varies by changing anisotropy 
of the hopping integral. 



V. COMPLETELY LOCALIZED 
EDGE/SURFACE STATES 

We have found that tight-binding models on the bipar- 
tite lattices, such as honeycomb and diamond lattices, 
have flat-band boundary states covering the whole BZ, 
when the anisotropy of their hopping integrals is suffi- 
ciently large. In general when systems have completely 
flat bands, one can construct a wave function which is 
spatially localized on a finite number of sites, as schemat- 
ically shown in Fig. [HI It is possible only when the fiat 
band covers the whole BZ. In this section, we show the 
localized wavefunction in the present models. This wave- 
function is exponentially decaying in the direction normal 
to the boundary, while on the outermost atomic layer, the 
wavefunction is nonzero only on a single site. We con- 
sider semi-infinite systems for the honeycomb lattice with 
the zigzag edge and those for the diamond lattice with 
the (111) surface, and similarly to the previous sections, 
the A sublattice faces the boundary. 

We first consider the honeycomb lattice in the half 
plane ^ < 0. The zigzag edge is along the x axis, and 
the origin O is set to be one site on the zigzag edge. The 
state at the site —ma.i — na.2 is denoted as \Amn) where 
m = (r?i, n) are integers, and O = (0, 0). We assume that 
for the fiat-band state the wavefunction at the boundary 
states is nonzero only at O. We express the wavefunction 
of the fiat-band states as 



I*) 



II «" 



(28) 



where amn is the amplitude for \Amn)' From the 
Schrodinger equation, we have the relations between am- 
plitudes as 



tia^i 



har, 



tsCtn 



0. 



(29) 



Intriguingly the solution of the above sequence is the 
same as the following problem. Consider a mover in the 
xy plane on the grid (Fig. [9|). The mover is first on the 
O = (0,0) site. At each step it moves by (1,0) with a 
probability Pi, by (0,1) with a probability P2 (Fig. [9|), 
and the movement is finished otherwise. Finally after 
N = m-\- n steps the probability Pmn that the mover is 



Surface 




FIG. 8: Schematic of spatial distribution of the flat-band 
surface states. The circles are the A sublattices, while the 
amplitudes on B sublattices are zero and are omitted. The 
dotted circles show that their amplitudes of the wavefunction 
are zero. The line thickness of the circle shows the magnitude 
of the amplitude. The top atom has the largest amplitude, 
and the distribution of the wavefunction spatially spreads to- 
ward the interior with exponential decay. This picture of the 
fully localized states are in common with the honeycomb lat- 
tice with a zigzag edge and the diamond lattice with the (111) 
surface. 



on (m, n) (m, n > 0) along the shortest path is given as 

(30) 
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mini 

By replacing the probabilities of the movement with the 
ratio of the hopping integrals, 



n^ 






(31) 



where i = 1, 2, we have the amplitude of the wavefunction 
as 



1 {m -\- n)l 
Zh mini 



■^) (m,n>0),(32) 



where Zh is the normalization constant. For m < or 
n < 0, amn vanishes. In this spatial representation of the 
wavefunction, the condition for the existence of the local- 
ized states on the boundary is that the wavefunction is 
normalizable. Generally, when ti^2,3(> 0) are independnt 
parameters, Zh satisfies the following relation: 
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(33) 



Therefore for ^^i^ < 1 the renormalization constant 
Zh converges; namely the completely flat-band states ap- 
pear. 

Next we consider the diamond lattice with the (111) 
surface in the half space x ^ y ^ z < 0. The surface are 
located onx + ?/ + z = plane, and the origin O is set 
to be (0, 0, 0). The wavefunction at mdi + nd2 + Zeds is 
denoted as \Arank)^ where d^ = r^+i — ti for i = 2,3,4. 
Then, the wavefunction of the flat-band state localized 
on the origin is denoted as 

1^) = Yl ^rnnk\Amnk). (34) 

m,n,/c 

where amnk is the amplitude for \Amn)- By the 
Schrodinger equation, we have the relations between am- 
plitudes as 
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Similarly to the case of the honyecomb lattice, the so- 
lution for the sequence is obtained by the shortest path 
problem in 3D. let Pi, P2 and P3 denote the probabil- 
ities for moving along (1,0,0), (0,1,0) and (0,0,1) re- 
spectively. The probability Pmnk of the shortest path 
problem from O to {m^n^k) is given as 



{m - 
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(36) 



By replacing the probabilities of the movement with the 
ratio of the hopping integrals, we have the amplitude 

^mnk as 



^mnk 



1 (tti + n + /c)! 
Zd mlnlkl 



(37) 



where Z^ is the normalization constant for the wavefunc- 
tion 1^), and this amplitude is nonzero only when tti, n, 
k are all nonnegative. The condition for the existence 
of the flat-band is that Z^ converges. Z^ satisfies as the 
following relation 
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(38) 



Therefore for ^2+^3+^4 <- i ^-j^^ renormalization constant 
Zd converges; the completly flat-band states appear. 

Those two convergences for Zh and Zd agree with 
the previous sections, where we use the representations 
of the Bloch wave functions. However, we note that 
Eqs. (|33])(|38|) give sufficient conditions for the conver- 
gence of the wavefunction. In Appendix [Aj we calculated 




FIG. 9: Schematics of the grids for the shortest path prob- 
lem. The mover is at first on the O, and it moves along e^ 
with the probability Pi, where (e^)^ = 6ij. 



models can be made in several ways: for example, split- 
ting the vertex of unipartite lattices. In addition when 
the system has anisotropy of the hopping integral along 
the direction normal to the surface, the boundary flat- 
bands appear. Thus to find candidate materials for 
the flat-band surface states, one need to find bipartite 
systems with strong anisotropy, and when the bulk is 
gapped, it may have the complete flat band depending 
on the surface orientation. 

The flat bands discussed in this paper have different 
properties compared with the conventional flat bandsi"— . 
Basically the conventional flat bands which cover the 
whole BZ originate from the difference of the numbers of 
the sublattices in the unit cell, while our system does not. 
Moreover, unlike those conventional flat-bands states, 
our flat-band states can be isolated from the bulk if the 
anisotropy becomes sufficiently large; in other words, the 
flat band is located in the middle of the bulk gap and the 
gap increases with the anisotropy. This property enables 
us to construct the fully localized boundary states. 



Zh_ and Zd ioT ti = t and the others being unity, where t 
is the hopping parameter of the anisotropy. The results 
agree with the conditions in Eqs. (|33]) (|38|) . From those 
results, we find that the edge/surface- localized wavefunc- 
tions spatially expand toward the interior of the systems, 
while the amplitude decays exponentially with the dis- 
tance from the surface (Fig. [8|). 



VI. SUMMARY 

Existence of flat-band boundary states is shown for 
tight-binding models on the diamond lattices. In the 
tight-binding model on the diamond-lattice model, we 
found that the distribution of the flat-band surface states 
at the Fermi energy is related with gapless loops in the 
bulk dispersion; the verge of the distribution of the flat 
bands in surface BZ is identical with the Fermi loop of 
the bulk projected on the surface direction. By increas- 
ing anisotropics of the hopping integrals, we showed that 
the distribution of the bulk Fermi surface shrinks and 
then disappears in the BZ. Then the surface flat band 
covers the surface BZ completely when the anisotropy 
is sufficiently large. In addition, we found a topologi- 
cal phase transition in the bulk Fermi loop in the tight- 
binding model on the diamond lattice by changing the 
anisotropy. These flat bands are understood topologi- 
cally within the theory by Ryu and Hatsugai^. Lastly, 
for strongly anisotropic cases we constructed wavefunc- 
tions which is localized only on one site on the boundary 
layer. From the Schrodinger equation, we showed that 
the spatially localized flat-bands states exists in the hon- 
eycomb and diamond lattice models. The wavefunction 
is localized at a single site in the outermost edge/surface 
layer, while it expands inward with exponential decay. 

In general, to obtain the localized boundary flat-band 
states, a bipartite lattice is necessary. Bipartite-lattice 
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Appendix A: Calculations for normalization 
constants 

In Appendix, we calculate the renormalization con- 
stants, Zh and Zd in Section |Vl for special cases where 
one of the hopping integrals is t and the others are unity. 
In the two cases, the hopping with t is perpendicular to 
the boundary. 



1. Honeycomb lattice 

For ti = t = x~^ and t2 = h = 1, Zh is given as 
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as long diS X = t~^ satisfies \x\ < ^. Therefore for |t| > 2 
there are flat-band states which are localized on one site 
on the edge in the honeycomb lattice with zigzag edge. 

2. Diamond lattice 
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and A^B G 3?. By changing the variable as tan/3 = v, 
the above equation becomes 
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where 6^ = 1 — 4- and c^ = 1 — \. The integral is further 
transformed to the elliptic integral of the first kind, 
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where t~^ = F for notational brevity and we assume 
< F <C 1. Later we can analytically continue the result 
with respect to Y to discuss the condition for convergence 
of Zd. i^ T] are complex numbers, and the path of the 
integrals are along the unit circle C : \z\ = 1. In the 
above equation, by using ^ = e*^ and rj = e^^ where 
^, (/) G 5R, we have 
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where the variable is changed as v 



bs 



By the 



above result, for < F < |, /c^ satisfies < /c^ < 1, and 
Zj converges. On the other hand, for F < 0, because Zj 

72 ._.__ f_ .1 < y < 0. 



is an even function of Y, Z^ converges for 



(A3) 



Namely, the condition for the convergence of the normal- 
ization constant Zd is \t\ < 3. 
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